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Abstract 

Two possible reasons for the narrow peak observation in the invariant mass spectrum of the pK^ close to 
1.54 GeV/c^ are discussed. 



1 Introduction 



Five different experiments have reported the observation of a narrow peak in the invariant mass spectra close 
to 1.54GeV/c^ with a width below 25MeV/c^. These peaks are considered as a signature of a five quark 
resonance, the so-called 6*+. 

Three photo-production experiments |21 El have observed the 0^ decay to 7iK+ in three- or four- body 
final state. The neutron was reconstructed from the missing mass and energy. 

The 9'^ decay to was observed in v^~{p^-) collisions with nuclei [J. This analysis has combined 
data samples from different bubble chamber experiments with a large range of neutrino momenta. 

Using a beam, the same channel of the 0^ decay was observed by the DIANA collaboration |S] . 
The Monte Carlo simulation for this experiment is rather simple for the following reason: first - the beam 
momentum spectrum is explicitly shown in the publication, second - all particles in the final state (p and 
K^^TT+TT^) were reconstructed and identified. 

The result of the simulation for the DIANA experiment are presented in the section below. One section 
further the result of the experiment "independent" simulation shows the possible source of the narrow peak 
in the invariant mass spectrum. 

2 Monte Carlo simulation. 

The experimental momentum spectrum of the incident beam particle K"'' (Fig.^ is shown for two cases: the 
momentum spectrum for all events with the identified (top) , and for events where both the K^? and the 
proton were identified [Hj (bottom). The momentum spectra look very different, both in momentum range 
and in shape. For the events with pK^ in the final state the momentum profile has an asymmetric "bell" -like 
shape with the average momentum close to 470MeV/ c. 

The K"*" momentum spectrum was simulated according to the experimental distribution shown in Fig.^. 
The momenta of the outgoing particles are generated uniformly in the available Lorentz-invariant phase 
space. Two processes are simulated : 

• three body final state - a nucleus play a role of the recoil particle: 
K+Xe ^pK^Xe' 

• two body final state - a charge-exchange reaction on the neutron of the Xe nucleus: 
K+n -^pKl 

The momentum is shared by two particles only. 

The pK^ invariant mass spectrum for three body final state simulation is shown in Fig.|23 . The compar- 
ison with the experimental distribution ("red" histogram) shows a reasonable agreement in positions, widths 
and shapes of both spectra. 

The invariant mass spectrum for the two body final state is shown in Fig.|2t. In contrast to the three 
body final state pK^Xe' the mass distribution resembles narrow peak with a sharp maximum at 1.55GeV/c^. 

The simple MC test with the K+ momentum spectrum of the gaussian shape (the mean at 450MeV/c 
and the width 15MeV/c) shows that the narrow momentum range of the incident beam is the reason 
why the narrow peak is seen in the mass spectrum in case of two body final state (Fig.|21;). The shift of 
the K+ momentum distribution by 20MeV/c gives the proportional shift by lOMeV/c^ in the mass peak 
position (Fig.|3i) and puts the maximum at 1.54GeV/c^. 

The simulation of the pK^ kinematic disagree with the observed invariant mass spectrum only in one 
point - peaks are shifted by 10-15MeV/c. This shift may come from the fact that in the MC the binding 
energy of the nucleus was not included in the simulation. There is an explicit indication in reference [S] 
that the beam momentum was shifted by ±15MeV/c in different data sample, but there is no indication 
about the corresponding statistics, accumulated with a particular beam momenta. As result, the MC events 
were not weighted and that may also cause some kind of shift in the peak position. The difference in the 
width is not too important -the MC shows that with statistic of about 30 event, the narrow peak may be 
observed even if the original distribution is wide (Fig.QJ. The narrow two-bin peak distribution, when at 
least 50% of entries belongs to two bins only, was observed in 2% cases - a non-negligible number. 



3 The kinematic "reflection" at 1.54GeV/c^. 



Since the mass region around 1.54GeV/c^has so much importance in the narrow peak search, another source 
of trouble must be pointed out. 

In a data sample of V'^ candidates {V'^ is used as a generic name for K^, A or A) the narrow peak in the 
mass spectrum may be created in the following way. 

A decays in two particles of opposite charge. Since the type of the particle is not known in the 
experiment the different mass hypothesis have to be assigned to each track to explore all combinations. 
With pion and proton/anti-proton masses assigned to corresponding tracks, each V'^ may represent K^, A 
or A. 

The most convenient way to present the kinematic of the decays is an Armentero-Podolansky plot [HI , 
shown in Fig.[31, where the three ellipses correspond to K°— >7r+7r~, A^p tt^ and A— i>p7r+decays. The ellipses 
have two overlap regions where this species are kinematically indistinguishable. Therefore, without an 
explicit cut on the corresponding mass hypothesis, there is always a cross contamination. For example 
may be not a "true" but a A or A in these narrow kinematic regions. 

In second step the V'^ is then combined with the same physical track which was already used to for the 
V'^ itself. Two class of combination may be defined: 

• True and the 7r+ of the but considered as a proton (class 1). 

• True A but considered as a K° and a positive track of the A which is a true proton (class 2). 

Figures El and [71 show the invariant mass distributions for the two classes of combinations defined above. 

The standard cut used in all experimental analysis is a cut on the invariant mass of the tt~^tt~ combination 
around mass (|m^o ±m^+^- | < 3-a). The usual mass resolution for two track combinations is limited to 
cr=2-4MeV/ c. For class 1 this cut does not change the invariant mass distribution (true K°), but for class 2 
this cut removes combinations with small masses and retains combinations around 1.54GeV/c^(Fig.|Hl). Fig.j^l 
shows the combined (class 1 and class 2) mass spectrum. 

It is hardly believable that a straight forward double track counting may happen in a real analysis. At 
least for DIANA this possibility is excluded by the nature of the experiment - in the bubble chamber all tracks 
are visible. The problem might be important for the electronic experiment if the cut on V'^ masses is not tight 
enough, thus the K° - A cross contamination is not removed completely, and the pool of reconstructed tracks 
is not free from "clones" or "ghosts" - de-facto the same track reconstructed twice due to , for example, some 
noise in the tracking detectors. The kinematic of such "clones" is very close to each other thus simulating 
the situation discussed above, naturally with a slight variation in the mass position and the width of the 
peak. 

4 Conclusion. 

The result of the simple MC shows that the invariant mass spectrum of the pK^combinations has a narrow 
structure close to 1.54GeV/c^ due to the narrow momentum range of the beam. If this is true, then further 
data analysis using the full moment range of the beam (see Fig.QJi) will cause the broadening of the observed 
peak and yield the answer about its nature. 

The artificial peaks in the invariant mass spectra due to kinematic reasons are well known. The only 
surprise is that such an artifact precisely coincides with the value of 1.54GeV/c^. 

The recent publication [T] shows that the nK~^ channel is not free from the kinematic reflections. 

Finally the question about the pentaquark signature will be found in a very classical way - increasing the 
statistic involved in the analysis and hopefully in a different experimental conditions. 
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Figure 1: a) The range and equivalent momentum (upper scale) of the K'*' beam for the events with 
identification only. The arrows indicate the momentum region of the p selection, b) The momentum of 
the beam for the events with p identification. The momentum distributions were taken from |S] 
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Figure 2: The experimental beam momentum [S] and MC mass spectra distribution corresponding to: b) 
reaction K+Xe — s-K^pXe'; c) reaction K+n ^K^p; d) the summ of both b) and c); The histogram in red 
corresponds to the experimental mass distribution from j^. 
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Figure 3: The beam momentum simulated according to the gaussian shape and MC mass spectra distribution 
corresponding to: b) reaction K+Xe — ^K^pXe'; c) reaction K+n — >K^p; d) the summ of both b) and c); The 
histogram in red corresponds to the experimental mass distribution from 0- 
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Figure 4: The invariant mass distribution for pK^ combination in case, when at least 50% of entries were 
observed in two bins (red) and the "true" distribution shape (black) in case of big number of entries. 
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Figure 5: The Armenteros-Podolansky plot for the V'^ candidates, pj^ and pt are the laboratory longitudinal 
and transverse momenta, respectively, of the decay tracks with respect to the y° direction. At the crossing 
point of the two ellipses the and A kinematically looks the same. 




Figure 7: The invariant mass distribution for pK^ combination of class 2. 
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Figure 8: The invariant mass distribution for pK^ combination of class 2 when the K^? mass is Umitcd within 
the Kg mass region. 
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Figure 9: The combined invariant mass distribution for pK^ combination of class 1 and class 2. when the 
Kg mass is limited within the K^ mass peak region. 



